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ABSTRACT 


INVESTIGATION  OF 
THE  PASCHEN  CURVE  OF  NITROGEN 
VIA  THE  APPLICATION  OF  NANOSECOND 
PULSED  ELECTROMAGNETIC  RADIATION 

David  W.  Scholfield 

B.S.,  Wichita  State  University,  Wichita  KS.,  1980 
M.S.,  University  of  New  Mexico,  1990 

In  order  to  obtain  lasing  in  microwave  excited  gases  a 
knowledge  of  the  Paschen  curve  is  required  so  that  estimates 
can  be  made  of  the  required  intensity  of  the  electromagnetic 
field  for  a  given  pulse  duration.  Data  in  the  region  below 


10'®  Torr-s 

with 

electromagnetic 

fields 

in 

excess  of 

lOkV/cm/Torr 

had 

been  unobtainable  due 

to 

technical 

limitations . 

If 

accurate  lasing 

research 

via 

microwave 

pumping  is  to  be  done  in  this  region,  experimental  data  must 
be  obtained  to  either  validate  or  repudiate  current  theory 
relating  to  this  region.  An  investigation  of  the  Paschen 
curve  of  nitrogen  from  pressures  of  0.336  Torr  to  685  Torr  via 
the  application  of  nanosecond  pulsed  electromagnetic  radiation 
has  been  investigated.  Breakdown  of  the  gas  was  achieved  via 
pumping  by  short  electromagnetic  pulses  and  verified  through 
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the  use  of  an  optical  multichannel  analyzer.  Previous 
investigations  have  utilized  monof requency .  electromagnetic 
signals.  Adjustments  for  the  short  pulse  nature  of  this 
investigation  have  been  derived  and  are  applied  to  the  Paschen 
curve.  Experiments  have  been  conducted  utilizing  a  pulser 
with  machine  parameters  of  1  J  delivered  during  1.375  ns  FWHM 
into  a  50  Q  coaxial  line,  at  a  peak  power  of  1  GW.  By 
investigating  this  previously  unexplored  region,  data  has  been 
acquired  in  the  region  of  the  relativistic  Paschen  curve  as 
hypothesized  by  Graham  and  Roussel-Dupre .  Experimental 
results  yield  no  observable  relativistic  effects  under  the 
conditions  of  this  experiment.  Additionally,  the  collisional 
frequency  limit  has  been  experimentally  identified. 
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I .  INTRODUCTION 


Experiments  at  the  Institute  for  Nuclear  Science,  Tomsk, 
Siberia,  in  the  former  USSR  have  demonstrated  an  interest  in 
the  ability  to  lase  diatomic  nitrogen  via  extremely  short 
microwave  pulses.  In  these  experiments,  lasing  was  achieved 
in  Nj  at  pressures  less  than  or  equal  to  one  atmosphere^- 
Excitation  of  the  gas  was  achieved  with  pulses  which  were 
approximately  between  16  and  50  ns  in  duration.  The  frequency 
of  the  microwaves  is  reported  to  have  been  3  GHz. 

To  lase  gas  with  microwaves  requires  a  knowledge  of  the 
Paschen  curve,  which  is  a  graphical  representation  of  the 
solution  of  the  gaseous  electrical  breakdown  equation. 
Breakdown  is  that  condition  in  which  the  final  electron 
density  due  to  some  form  of  excitation  is  10®  times  the 
initial  electron  density.  Knowledge  of  the  Paschen  curve 
allows  estimates  to  be  made  of  the  required  intensity  of  the 
electromagnetic  field  for  a  given  pulse  duration.  Data  in  the 
region  below  10'®  Torr-s  with  electromagnetic  fields  in  excess 
of  lOkV/cm/Torr  had  been  unobtainable  due  to  technical 
limitations.  If  accurate  lasing  research  via  microwave 
pumping  is  to  be  done  in  this  region,  experimental  data  must 
be  obtained  to  either  validate  or  repudiate  current  theory 
relating  to  this  region.  This  document  describes  an 
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experiment  preformed  in  this  region. 


Figure  1,  see  page  4,  is  a  plot  of  the  various  Paschen  curves 
of  previous  researchers  available  in  the  literature.  The 
abscissa  represents  PXp,  which  is  the  pulse  duration 
multiplied  by  the  operational  pressure,  while  the  ordinate 
represents  E^/P,  which  is  the  applied  electromagnetic  field 
divided  by  the  operational  pressure. 

The  first  two  lines  on  the  upper  right  hand  side  of  figure  1 
are  predications  based  on  the  theory  of  Ali  and  Coffey^  for 
monofrequency  microwave  excitations  of  the  upper  atmosphere  at 
10  GHz  and  3  GHz  respectively.  For  this  study  Ali  and  Coffey 
assumed  that  the  species  from  which  ionization  occurred  was 
oxygen,  which  has  a  first  ionization  energy  of  12.5  eV. 

The  third  line  from  the  top  of  the  upper  right  hand  side 
predicts  that  level  above  which  breakdown  must  occur  based  on 
data  acquired  by  Sullivan  et  al.^  with  monofrequency 
microwaves  excitation  of  air  at  9.6  GHz.  Sullivan  utilized  a 
10  to  200  kW/cm^  gyrotron  with  pulses  of  between  3  and  30  ns. 

The  third  line  from  the  bottom  left  hand  side  is  a  prediction 
based  on  the  theory  of  MacDonald^  for  monofrequency  microwaves 
excitation  air  at  3  GHz. 
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The  second  line  from  the  bottom  left  hand  side  is  based  on 
data  acquired  by  Felsenthal  and  Proud^  for  a  pulsed  dc 
experimental  system  utilizinq  nitrogen  as  the  working  medium. 
The  pulsed  system  was  capable  of  generating  a  4  to  30  kV  pulse 
with  rise  times  a  short  as  0.3  ns. 

Finally,  the  bottom  line  on  the  left  hand  side  is  based  on 
data  acquired  by  Tetenbaum  et.  al.®  which  measured  breakdown 
due  to  a  pulsed  S-band  microwaves  with  pulse  durations  of 
between  100  to  600  ns.  The  pressure  at  which  Tetenbaum 
operated  varied  from  1  to  1000  Torr. 

The  variation  in  the  placement  of  the  theories  and  the  data 
over  approximately  two  decades  is  due  to  the  different  media 
and  different  frequencies  of  investigations.  Generally, 
higher  frequencies  reside  closer  to  the  upper  right  hand 
corner  of  figure  1.  This  is  due  to  the  effect  of  frequency  on 
the  calculation  of  the  effective  electric  field.  Higher 
frequencies  tend  to  transfer  less  energy  to  an  electron. 
Likewise,  nitrogen  gas  tends  to  reside  closer  to  the  bottom 
left  hand  side  of  figure  1  with  air  and  oxygen  residing  closer 
to  the  upper  right  hand  side. 


3 


Compilation  of  Previous  Research 


Figure  2  on  page  7  is  a  plot  of  pressure  versus  pulse  length 
compiled  from  a  survey^'^^  of  other  research  -  R.F.,  spark  gaps, 
etc.  -  relevant  to  this  problem.  The  darkened  areas  are  those 
which  have  been  investigated,  the  boxed  region  is  the  area  of 
interest  to  this  investigation. 

Utilizing  a  recently  constructed  pulse  power  device  designated 
H-Two,  it  was  possible  to  produce  the  required  field  strengths 
and  pulse  lengths  necessary  to  perform  this  experiment.  The 
H-Two  device  is  capable  of  producing  a  1.375  ns  FWHM  pulse  of 
up  to  720  kV  across  an  open  transmission  line.  Breakdown  of 
nitrogen  was  achieved  via  excitation  by  broadband 
electromagnetic  pulses,  with  a  peak  frequency  content  at 
approximately  200  MHz. 

This  investigation  has  examined  the  Paschen  curve  of  nitrogen 
from  pressures  of  0.336  Torr  to  685  Torr  via  the  application 
of  nanosecond  pulsed  electromagnetic  radiation.  This  extends 
the  efforts  of  earlier  researchers  by  investigating  the  region 
in  excess  of  10  kV/cm/Torr  for  E^/P  and  less  than  10'®  Torr-s 
for  PXp,  which  had  been  experimentally  denied  due  to  technical 
limitations.  An  extensive  literature  search  was  unable  to 
uncover  any  experimental  data  in  this  region. 

The  previous  efforts  utilized  monofrequency  sources,  whereas 
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this  effort  utilized  a  broad  band  source  with  a  frequency 
content  peak  at  approximately  200  MHz.  Since  the  results  of 
Felsenthal  and  Proud  were  obtained  for  a  pulsed  D.C. 
experimental  system,  and  the  results  of  Ali  and  Coffey  were 
calculated  for  10  GHz  and  3  GHz,  the  broadband  signal  to  be 
bounded  on  the  bottom  by  the  work  of  Felsenthal  and  Proud,  and 
on  the  top  by  the  results  of  Ali  and  Coffey. 
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Figure  2.  Pressure  vs.  Pulse  Length  Compilation 


Investigation  of  the  region  in.  excess  of  10  kV/cm/Torr  for 
Eo/P  and  less  than  10'®  Torr-s  for  PXp  has  yielded  data  in  the 
regime  of  the  relativistic  Paschen  Curve  as  theorized  by 
Graham  and  Roussel-Dupre^®,  see  figure  3  on  page  9.  Graham  and 
Roussel-Dupre  developed  the  relativistic  theory  to  describe 
phenomenon  due  to  lightning  in  the  upper  atmosphere.  The 
theory  was  developed  assuming  intense  electromagnetic  fields, 
with  electron  flow  both  perpendicular  and  parallel  to  the 
propagation  of  the  electric  field. 

This  document  is  divided  into  four  sections.  Section  I,  just 
presented,  was  the  introduction,  which  discussed  the 
motivation  for  this  research,  the  previous  investigations  of 
prior  researchers,  and  the  uniqueness  of  this  experimental 
investigation.  Section  II  develops  classical  and  relativistic 
theories  for  broadband  electromagnetic  radiation,  extending 
previous  classical  theories  for  monofrequency  electromagnetic 
sources .  Section  III  describes  the  pulse  power  device  used  to 
generate  the  broadband  electromagnetic  pulse,  the  test 
chamber,  and  the  diagnostics  used  to  monitor  the  pulse. 
Section  III  continues  with  a  discussion  of  the  procedure  to 
acquire  the  data,  the  process  used  to  reduce  the  data,  and  the 
analysis  perfo2rmed  on  the  data.  Finally,  Section  IV  presents 
the  conclusions  derivable  from  the  analysis  performed  on  the 
data . 
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II  A.  CLASSICAL  THEORY 


This  section  develops  classical  theory  ,  for  broadband 
electromagnetic  radiation.  This  development  represents  an 
extension  to  current  theories,  which  have  been  developed 
predominately  for  monofrequency  electromagnetic  excitation. 

For  nanosecond  pulsed  electromagnetic  radiation,  minor 
adjustments  to  monofrequency  classical  theory  must  be  made. 
As  has  been  shown  by  previous  authors^'^'^^,  the  equation  which 
governs  the  time  rate  of  growth  of  breakdown  is  the  continuity 
equation, 

{Dn)  (D 

at 

where:  n  =  electron  density 

Vi=  ionization  rate 
Va=  attachment  rate 

D  =  electron  diffusion  coefficient 
When  an  electromagnetic  pulse  is  short  compared  to  the 
diffusion  time,  the  electron  diffusion  coefficient  may  be 
neglected®.  This  allows  equation  1  to  be  rewritten  as 

an-Avu  (2) 


where :  Av  =  Vj  -  v^ 


10 


Solving  equation  2,  assuming  Av  is  not  a  function  of  time, 
yields 

In— “AvTo 


where:  n^  =  final  electron  density 

n^  =  initial  electron  density 
Xp=:  pulse  length 

Gould  and  Roberts^^  established  the  generally  recognized 
breakdown  ratio  of  nj/n^  as  10®. 

lnlO®-l8.4“AvTp 

For  breakdown  to  occur  it  is  required  that  Vi»Va.  Thus 

^i-  18.4  (5) 

P  P^p 


where:  P  =  gas  pressure 

Following  the  development  of  Lupan^®,  the  functional  dependence 
of  the  ratio  of  vjV  to  Eeff/P  can  be  derived. 
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The  average  value  of  this  ratio. is 


where : 


Assuming 


where^® 


P  P 


2E 


(6) 


N  =  gas  density 

a(E)  =  collisional  cross  section  of  the  gas  molecule 

E  =  energy 

lUg  =  electron  mass 

<CT(E)  *  [2E/me]  =  average  collisional  rate 

coefficient  integrated  over  a  distribution  function, 
a  Maxwellian  distribution  for  E, 


2E  2 


E 


exp{  „ 
3 


-3E\  dE 


[|e] 


(7) 


e- 


—ksT+ 
2  ® 


and  kg  =  Boltzmann's  constant, 

T  =  temperature, 

Q  =  electron  charge, 

Eeff  =  effective  E-field, 

5eff  =  effective  fraction  of  energy  transferred 
from  an  electron  to  a  molecule  in  a 
collision, 

Vgff  =  effective  collision  frequency. 
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Equation  7  can  be  rewritten  as 


2 

,/jrN 


2e  2eJ  P 

El 


— /Fexp{ 


-3E 

2e 


(8) 


here  Ei  is  the  first  ionization  potential  of  nitrogen.  The 
functional  dependence  of  the  first  three  terms  contained  in 
the  integral  can  be  estimated^®  for  the  field  values  of  this 
investigation  based  on  data  of  the  ionization  cross  section  of 
nitrogen,  and  assuming  that  ionization  is  a  single  event 
process.  For  E  <  400  eV  this  estimate  yields, 

— o(£?).^^-6xI0’a^/F(£-Pi)exp(^^^)  (9) 
P  \  PjK 


where^’:  a=0.2  cm‘^  Torr'^ 

P,2=142.5  eV 
Ei=15.5  eV 


For  E  >  400  eV,  the  functional  dependence  of  the  first  three 
terms  contained  in  the  integral  are  estimated  to  be 


Zo(P).  (10) 

P  N  ^ 
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Substituting  equations  9  and  10. into  equation  8  yields, 


400 


E,-E, 


<Zi>-6xlO\|  f  e(E-E^)  exp{-^^}exp{^ 

P  ^271  J  Pn2 


}dE4 


10' 


J/ffexp  ( 


400 


-3E 

2e 


)d£] 


(11) 


Performing  the  integration  (refer  to  appendix  A)  and  combining 
with  equation  5  yields 


18.4 

Pi. 


exio 


\  2it 


27  .-3/2{ 


2«Pjgg 

2e +3 


8(PfBg)^ 

{26+3P|J2^  ^ 


2B+3p^ 


+e|]  exp  ( 


-3Ej 

2e 


) 


2«P«zggi 

2e +3  p^ 


2Pi,zg 

2b+3  Pj^ 


+Ei)exp( 


-3^1 

2b 


) 


2tce 

3 


exp 


zi2£)  (^+800)} 

B  3 


(12) 
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Which  reduces  to 


18.4 

Plr. 


6x10’’ 


e"^/^exp  ( 
2tz 


^  I  16g  ( pj^igE )  ^  ^  4tt  (Pj^gE ) 

26  ^  (2e  +  3P„2)3  (2e  +  3Pj,2)2 


+  10' 


N 


2j!lexp(3Sd££E.)  (^*800)) 
3  2e  3 


(13) 


For  e«pN2/  the  above  solution  approximately  yields 


Pt« 


2ti  2e  27  9 


+10' 


\  3 


2^16  -_/  -1200  +  3.Fi  2e 

exp  ( — 


2e 


)  (.^+800)} 


(14) 


For  e=PN2/  equation  13  yields 


Pv, 


\  271 


26  125 
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+10= 


N 


271E _ ,  -1200+3^i  2e 

— —  exp  ( - - )  ( 


2e 


h800)} 


(15) 
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For  e»pN2  and  Ei,  equation  13  approximately  yields 


—  e'^^^exp(— ^^){2apL+«PL^i 
PZp  N  2%  2e 


2ne 

3 


exp( 


-1200+3^i 

2e 


)  (^  +  800)} 


(16) 


Recall  that, 

e-ljc^rt  (17) 

This  discussion  will  assume®  that  Veff=5.3xl0  P  where  Vg^f  has 
units  of  s-^  and  P  has  units  of  Torr,  and  from  Appendix  B  that 
Ee££=0.310Eg  for  the  area  of  interest  for  this  investigation. 
Observing  that  the  average  electron  energy  in  a  microwave 
discharge  is  much  greater  than  (3/2)kBT,  that  mg=  9.11x10  gm, 
and  that^®  5gff=l .  3x10'^  equation  17  yields 

(0.310)^(1.602x10'^^)  10“^  (^^0)2 

1.3x10-2(5.3x10^)2(9.11x10-28)  P 


e-4. 63x10-^  ( 


OE, 


(18) 
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Equation  13  is  plotted  in  figure  4  on  page  18  with  the 
substitution  of  equation  18,  and  the  values:  a=0.2  cm  ^  Torr 
P„2=142.5  eV,  and  Ei=15.5  eV.  The  ordinate  represents  EJP 
while  the  abscissa  represents  PXp.  The  graph  of  this  Paschen 
curve  demonstrates  a  departure  from  the  traditional 
monofrequency  curves  of  previous  researchers .  At 

approximately  1x10'^  Torr-s  and  1000  V/cm/Torr  the  Paschen 
curve  reverses  and  transitions  to  the  collisional  frequency 
limit .  The  features  of  reversal  and  transition  to  the 
collisional  frequency  limit  have  not  been  addressed  in  the 
Paschen  curve  investigations  of  previous  researchers. 


II  B.  RELATIVISTIC  THEORY 


This  section  describes,  in  brief,  the  relativistic  theory  of 
G2rahain  and  Roussel  “Dupre  and  then  develops  a  one  dimensional 
theory  for  intense  electric  fields,  utilizing  classical 
scattering  theory  with  relativistic  corrections. 

Scattering  at  classical  energies  allows  the  electric  fields  of 
the  incoming  particles  to  be  treated  as  isotropic  fields  (see 


electron  is  capable  of  obtaining  a  velocity  sufficient  to 
exhibit  relativistic  effects^®.  When  this  occurs,  the  electric 
field  of  the  incoming  relativistic  particle  redistributes 
itself,  as  seen  by  the  scattering  center  (fixed  frame) . 
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The  equation  describing  this  redistribution  is 


_ 0 _ 

(1-Preisin^<l>) 


Where:  Q  =  Electron  charge  (1-6x10'^®  C) 

R  =  charge  separation 
Prel  =  v/c 

y  =[  1  - 

The  field  is  in  the  radial  direction,  but  the  electric  lines 
of  force  are  distributed  isotropically  only  when  Prei*0  •  For 
Pj.gi<#>0,  the  electric  force  is  reduced  by  a  factor  of  when 
viewed  along  the  axis  of  motion,  while  when  viewed 
perpendicular  to  the  axis  of  motion  the  electric  field  is 
increased  by  a  factor  of  y.  This  effect  is  shown  in  figure  6 
for  y~l  and  y~3  on  page  21,  with  the  magnitude  of  the  electric 
field  represented  by  line  length.  Under  this  influence  the 
effective  collision  frequency  has  been  shown®^  to  be 


rel 


Noting  equation  B-1  of  appendix  B  it  is  seen  that  the 
effective  E-field  is  a  function  of  the  effective  collision 
frequency. 


Eeff- 


[2(1  + 


Ci>^ 


)] 


1/2 
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Figure  6.  Electric  1 
Y=1  anc 


2: 


Thus  the  relativistic  ratio  Prei  is  a  function  of  Egff/P,  which 
itself  is  a  function  of  which  in  turn  is  a  function  of 

Thus  an  existential /transcendental  set  of  equations  must 
be  solved.  This  problem  has  been  addressed  via  a  numerical 
approach  by  Graham  and  Roussel-Dupre^^  for  air.  Figure  7,  on 
page  23,  is  a  plot  of  the  Paschen  curve  derived  by  Graham  and 
Roussel-Dupre  based  on  their  numerical  approach.  The  ordinate 
represents  E^/P  in  V/cm/Torr  while  the  abscissa  represents  PXp 
in  Torr-s.  The  solid  line  represents  the  relativistic  Paschen 
curve  for  air.  Features  of  interest  include  the  reversal  of 
the  curve  from  IxlO'^  V/cm/Torr  to  2x10''  V/cm/Torr,  and  the 
reestablishment  of  the  curve  above  3x10^  V/cm/Torr.  The 
hypothesis  put  forth  by  Graham  and  Roussel-Dupre  as  to  the 
cause  of  this  reversal  is  that  as  an  electron  achieves  higher 
velocities  the  electric  field  of  the  electron  becomes 
redistributed,  as  viewed  in  the  laboratory  frame,  due  to 
Lorentz  contraction.  As  the  electron  achieves  velocities 
approximately  equal  to  the  velocity  of  light  the 
redistribution  becomes  essentially  fixed,  that  is  to  say  that 
additional  energy  added  to  the  electron  can  no  longer  provide 
as  great  an  acceleration  as  before,  and  the  Paschen  curve 
reestablishes  the  original  behavior. 
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The  Paschen  curve  derived  by  Graham  and  Roussel-Dupre  assumed 
intense  E-fields,  B-fields,  and  electron  flow  in  two 
dimensions.  The  experimental  configuration  -of  this 
investigation  only  allowed  intense  E-fields  with  electron  flow 
in  one  dimension.  This  necessitated  extension  of  the 
classical  theory  into  the  relativistic  domain.  Appendix  C 
develops  a  relativistic  derivation  of  the  Paschen  curve 
beginning  from  first  principles.  The  results  of  this 
derivation  are  presented  in  figure  8  on  page  25.  The  ordinate 
represents  E^/P  in  V/cm/Torr  while  the  abscissa  represents  PXp 
in  Torr-s.  The  solid  line  represents  the  relativistic  Paschen 
curve  for  air,  the  long  dashed  line  represents  the  classical 
curve  for  nitrogen,  the  short  dashed  line  represents  the 
relativistic  curve  for  nitrogen  as  derived  in  Appendix  C. 
Note  that  all  three  theories  produce  a  reversal  in  the  paschen 
curve.  In  the  classically  derived  theory  the  reversal 
transitions  into  a  vertical  barrier.  In  the  relativistic 
theory  of  Graham  and  Roussel-Dupre  the  reversal  transitions 
through  yet  another  reversal  to  reestablish  the  Paschen  curve 
at  a  higher  value  of  E^/P.  The  relativistic  theory  developed 
in  Appendix  C  produces  a  reversal  which  continues  into  the 
upper  right  hand  quadrant  of  the  plot,  which  is  attributed  to 
the  one  dimensional  nature  of  the  theory  (  Refer  to  Appendix 
C  for  details  ) . 
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The  H-Two  device,  to  be  discussed  in  the  Experimental 
Configuration  section,  has  produced  electromagnetic  pulses  as 
part  of  this  investigation  capable  of  exploring  the  region 
below  10'®  Torr-s  with  electromagnetic  fields  in  excess  of 
lOkV/cm/Torr .  The  results  of  this  investigation  follow. 
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Ill  A.  EXPERIMENTAL  CONFIGURATION 


Section  III  describes  the  H-Two  pulse  power  device  used  to 
generate  the  broadband  electromagnetic  pulse,  the  test  chamber 
in  which  the  discharge  occurred,  and  the  diagnostics  used  to 
monitor  the  event.  This  section  continues  with  a  discussion 
of  the  data  acquisition  procedure,  the  process  used  to  reduce 
the  data,  and  the  analysis  performed  on  the  data. 

The  H-Two  pulse  power  device  was  used  to  provide  the  transient 
electromagnetic  events .  A  schematic  of  this  device  is  shown 
in  figure  9 . 


This  device  is  capable  of  developing  a  1  J,  1.375  ns  FWHM 
pulse  of  up  to  360  kV  into  a  50  Q  load,  or  720  kV  into  an  open 
circuit.  A  typical  oscilloscope  trace  is  shown  in  figure  10 
on  page  28.  The  abscissa  is  time  in  nanoseconds  while  the 
ordinate  is  volts.  Visible  is  the  rising  edge  of  the  pulse, 
the  occurrence  of  breakdown,  the  ringing  of  the  circuit  after 
the  arc,  and  the  exponential  decay  envelope. 
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Figure  11  on  page  3  0  is  a  plot  of  a  hypothetical  1  ns  duration 
triangular  pulse  with  a  peak  amplitude  of  1  MV  at  a  pressure 
of  10  Torr.  The  abscissa  is  Torr-s  while  the  ordinate  is 
V/cm/Torr.  The  solid  line  represents  the  relativistic  Paschen 
curve  while  the  dashed  line  represents  the  track  of  the  pulse. 
This  plot  demonstrates  that  a  machine  with  H-two's  performance 
can  investigate  the  nose  region  of  the  relativistic  Paschen 
curve . 

The  high  voltage  power  supply  was  a  Hippotronics  Labtrol 
capable  of  25  kV  charge.  Triggering  was  provided  by  a  Maxwell 
Laboratories,  Inc.  40127-B  Trigger  Channel,  which  in  turn 
drove  a  Northstar  Research  Corp.  Thyratron  Driver,  and  an 
Optical  Multichannel  Analyzer.  The  thyratron  driver  provided 
the  trigger  to  the  thyratron  switch  between  the  capacitors  and 
the  primary  inductor  as  shown  in  the  schematic. 

A  transparent  test  chamber  capable  of  pressures  from  100  mTorr 
to  7600  Torr  (gauge)  was  constructed,  and  clamped  directly  to 
the  H-Two  device.  The  chamber  was  a  76  cm  long,  30  cm  I.D. 
plexiglass  tube  with  a  thickness  of  2.54  cm.  The  end  through 
which  the  transmission  line  entered  the  test  chamber  was 
sealed  with  a  2.54  cm  thick  aluminum  bulkhead.  The  end  at 
which  the  parallel  plate  transmission  line  terminated  was 
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Figure  11 .  Triangular  Pulse  Track 


sealed  with  a  4.78  cm  thick  plexiglass  bulkhead  and  a  5.08  cm 
plexiglass  access  panel.  See  figure  12. 


Figure  12.  Test  Chamber. 

The  chamber  contained  a  section  transitioning  from  the  50  Q 
coaxial  transmission  line  to  the  parallel  plate  transmission 
line.  The  lower  plate  of  the  parallel  plate  transmission  line 
was  the  electrical  ground.  Two  center  conductors  were 
machined  for  the  experiment.  A  7.23  cm  wide  center  conductor 
was  manufactured  to  produce  a  parallel  plate  transmission  line 
with  an  impedance  of  50  £2  at  a  1  cm  AK  gap.  A  1  cm  wide 
tapered  center  conductor  was  manufactured  to  produce  a 
parallel  plate  transmission  line  with  an  impedance  of  50  Q  at 
a  0.14  cm  AK  gap. 
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The  data  acquisition  system  (see  figure  13  on  page  33) 
consisted  of  a  Zenith-248  286  computer  equipped  with  a 
National  Instruments,  Inc.  PC2A  GPIB  board.  The  software  used 
to  operate  the  data  acquisition  system  was  DCSOl,  version 
2.05,  by  Tektronix,  Inc.  The  286  computer  was  connected  via 
the  GPIB  board  to  a  Tektronix,  Inc  MUX  16  Multiplier.  The  mux 
box  was  required  so  that  the  software  could  communicate  with 
three  DCSOl  cameras.  The  cameras  were  mounted  on  the  front  of 
three  Tektronix  7103  oscilloscopes.  The  7103  oscilloscopes 
were  equipped  with  7A29  preamps,  and  7B10  time  bases  to 
provide  1  GHz  bandwidth  and  0.35  ns  risetime  capability  for 
signal  diagnostics. 

Optical  diagnostics  were  provided  by  an  E.G.&G.  Princeton 
Applied  Research  Model  1460  Optical  Multichannel  Analyzer  III. 
The  OMA  was  used  to  verify  that  a  N2  spectrum  was  being 
produced  by  the  discharge.  The  optical  data  was  acquired  via 
two  lengths  of  ultra  violet  rated  Mitsubishi  ST-U200D-FV  fiber 
optic  cable.  A  25  m  long  cable  viewed  the  discharge 
perpendicular  to  the  arc  while  a  50  m  long  cable  viewed  the 
discharge  down  the  axis  of  the  arc. 
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Figure  13.  DAS  Block  Diagram. 


The  data  lines  consisted  of  five  RG-214  cables,  with  N-type 
connectors,  running  from  the  screen  room  to  the  experimental 
area.  Inside  the  screen  room  the  cables  transitioned  to  152.4 
cm  long  RG-141  cables  with  SMA  connectors.  At  the 
experimental  area,  additional  lengths  of  RG-214  were  attached 
to  the  cables  to  provide  a  transition  from  the  experimental 
area  patch  panel  to  the  experimental  device.  All  data  lines 
were  then  measured  for  attenuation  characteristics  and 
electrical  length.  Since  the  coaxial  lines  were  operating  in 
TEM  mode  dispersion  was  negligible.  The  three  data  lines  most 
closely  matched  in  electrical  length,  and  with  the  best  match 
to  the  theoretical  attenuation  curve,  were  selected  as  the 
primary  data  lines  for  the  pulse  power  diagnostic, 
transmission  line  sensor,  and  the  tip  sensor.  Of  the  two 
remaining  data  lines,  one  was  utilized  as  a  trigger  line  for 
the  thyratron  drive  while  the  other  was  left  idle  to  provide 
a  backup  in  the  event  of  a  failure. 

Pressure  inside  the  test  chamber  was  sensed  by  two  Granville- 
Phillips  275  Convection  Gauges  and  by  one  Varian  865  gauge 
equipped  with  a  053  6  vacuum  transducer.  This  duality  provided 
for  a  cross-check  to  ensure  that  the  gauges  were  functioning 
properly. 
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The  pulse  power  switch  sensor  consisted  of  copper  clad  RG-141 
with  the  center  conductor  extending  0.165  cm  radially  inward 
from  the  outside  wall  of  the  gaseous  switch.  See  figure  14. 


rG-141  Center  Conductor 


Figure  14.  Pulse  Power  Switch  Sensor. 

The  outside  wall  represented  electrical  ground  and  was  in 
electrical  contact  with  the  copper  outer  conductor  of  the 
probe.  This  allowed  the  sensor  to  function  as  a  half  wave 
dipole  antenna  at  a  frequency  of  approximately  45  GHz. 
Equipment  to  provide  a  rigorous  calibration  of  the  pulse  power 
switch  sensor  was  not  available.  Therefore  an  estimate  of  the 
sensor's  signal  coupling  ability  was  calculated  based  on  the 
geometry  of  the  gaseous  switch.  The  gap  from  the  switch 
channel  to  the  grounded  wall  containing  the  sensor  was  8.839 
cm  in  length. 
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Attenuation  due  to  geometrical  considerations  is 

Attenuation  due  to  a  mismatch  between  the  transmission  line 
impedance  and  the  pulse  power  switch  sensor  radiation 
resistance,  Ra,  must  also  be  accounted  for.  The  radiation 
resistance  for  an  infinitesimal  dipole  antenna  has  been 
calculated  to  be^® 

i?^-80ir2(^)2 

which  for  the  geometries  and  wavelengths  under  consideration 
yields 

i2^-807l2(Jl^^)2-o.00096Q 


Since  the  transmission  line  and  the  data  cable  are  both  50  SI 
it  is  necessary  to  account  for  reflections  at  both  interfaces. 
The  voltage  reflection  coefficient,  F,  is 

0  •  00096~50 _ Q  99955 

0.00096+50 

This  yields  a  functional  attenuation  of  0  dB  at  each 
interface.  The  total  estimate  for  the  attenuation  of  the 
pulse  power  switch  sensor  is  34.5  dB. 
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The  transmission  line  sensor  consisted  of  a  capacitative  probe 
of  approximately  2200  pF,  constructed  with  a  2.54  cm  diameter 
disk  of  silver  coated  0.0005  cm  thick  mylar  (see  figure  15 
below) .  The  probe  was  recessed  0.737  cm  into  the  wall  of  the 
50  Q.  coaxial  transmission  line.  Equipment  to  provide  a 
rigorous  calibration  of  the  transmission  line  sensor  was  not 
available.  As  this  sensor  was  used  as  a  qualitative  health 
check  of  the  pulse  power  device,  a  rigorous  attenuation 
estimate  was  not  pursued. 


Silver  Coated 
0.0005  Mylar  Film 


Figure  15.  Transmission  Line  Sensor. 
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•  of  an  SMA  bulkhead  feedthrough  jack 

The  tip  sensor  consisted  of  an  bNA  ^  ^ 

^  fi.er  optic  poPt.  io.  the 
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(see  figure  16)  . 


Figure  16.  Tip  Sensor. 

;:,f  the  tip  of  the  ground  plane  of 

The  tip  sensor  was  mounted 

The  sensor  was  mounted  such 
the  50  Q  flat  plate  transttxon.  The 

ground  plane. 

.he  asse^ly  with  the  7.23  c.  wiae  oenten  conauotor  and  1  - 
„,ae  oenten  conaucton  was  then  calihnatea  fot  oonpUn.  a 

Phase  effects  with  an  HP  S510  vector  network  analyz  ^ 

oofficient  had  an  attenuation  of 
overall,  the  coupling  coef 

,  f-ioiir-Ps  17  and  18  on  pages  39  ana 
approximately  -25  dB  (see  figures 

respectively) . 
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Magnitude 


*-  •  r^haracteristics  of  the  fiber  optic  cables  are 
The  attenuation  characteristiu 

^  AO  a  nlot  of  attenuation  versus 

shown  in  figure  19  on  page  42,  a  plot  o 

•  r.  T  ra-irfal  pnath  in  nanometers  while 
wavelength.  The  abscissa  is  wavelength 

the  ordinate  is  attenuation  in  decibels. 

The  spectrograph  of  the  OMA  was  a  darrell-Ash  MonoSpec  27 
Spectrograph  with  an  entrance  slit  of  25  \i,  and  a  grati  g 
150  g/.™  blaze.  Ibe  output  of  the  Spectrograph  was  detecte 

by  an  E.G.&G.  Model  1254  SIT  CCD. 

sensitivity  curve  of  the  CCD  is  shown  in  figure  20  on  page 
43,  a  Plot  of  sensitivity  versus  wavelength.  The  abscissa  rs 
wavelength  in  nanometers  while  the  ordinate  is  sensit 
counts  per  photon. 

.he  CCD  was  controlled  by  an  E.C..G.  Model  1216  Multichannel 
Detector  Controller  and  was  triggered  by  and  E.G..G.  Model 
1211  High  Voltage  Pulse  Generator. 
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ation 


.S,G.  Model  1254  SIT  CCD 


Ill  B.  EXPERIMENTAL  PROCEDURE 


The  experimental  procedure  utilized  the  pulser's  natural  pulse 
length  of  1.375  ns  FWHM.  The  charge  voltage  was  set  at  18  kV 
Yielding  a  peak  signal  of  approximately  320  kV  on  the  50  Q 
coaxial  transmission  line  leading  to  the  test  chamber.  Upon 
reflection  of  the  pulse  at  the  open  circuit  termination  of  the 
50  a  parallel  plate  transmission  line,  the  voltage  of  the 
pulse  doubles  to  640  kV  across  a  1  cm  gap.  No  seed  ionization 
was  employed  in  this  experiment . 

Table  1  shows  the  pressure,  number  of  shots,  and  number  of 
shots  yielded  by  the  experimental  procedure.  The  annotation 
in  the  table  of  an  evacuation  indicates  that  the  test  chamber 
was  evacuated  to  approximately  100  mTorr  and  then  backfilled 
with  dry  nitrogen  to  the  desired  pressure.  The  test  chamber 
was  reconfigured  with  the  narrow  center  conductor  and  the  AK 
gap  adjusted  to  0.14  cm  to  provide  a  50  £2  parallel  plate 
transmission  line  for  the  0.336  Torr  shots.  The  number  of 
shots  yielded  indicates  the  number  of  shots  in  which  an  arc 
has  been  detected.  The  procedure  for  detecting  an  arc  is 
discussed  in  the  next  section.  Data  Reduction. 
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Evacuated 

0.930 

35 

11 

18* 

Evacuated 

685 

295 

150 

73 

29.5 

19.5 

1.07 

12 

14 

10 

12 

11 

22 

13 

10 

10 

10 

10 

10 

7 

5 

C 

18 

18 

18 

18 

18 

18 

18 

18 

0.506 

19 

D 

20.125 

0.506 

21 

f 

Evacuated 

15 

Q 

18 

2.08 

C 

18 

1.22 

11 

D 

c 

18 

1.01 

10 

D 

a 

18 

0.890 

10 

b 

18 

0.790 

12 

i 

•1  A 

18 

0.700 

13 

10 

o 

18 

0.590 

12 

o 

Evacuated 

3.4 

4.8 

7.0 

33 

23 

33 

18 

5 

7 

18 

18 

18 

Evacuated 

0.336 

44 

24 

20 

Table  1 .  pressure  versus  Nuinber  o£  Shots 
685  Torr  to  0.336  Torr 


V.  this  shot  sequence 

*  Transformer  rebuild  after  tni 
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IXX  C .  DATA  REDUCTION 


The  data  was  ac,uired  as  described  in  the  experimental 
procedures  section.  The  data  was  stored  by  the  DCSOl  software 
in  log  files  which  allowed  the  DCSOl  software  to  marntarn 
large  amounts  of  data  on  files  of  no  greater  than  2200  bytes. 
These  files  had  to  be  converted  to  ASCII  files  by  a  DCSOl 
software  utility  before  data  analysis  could  begin.  After 

^  acipTi  the  data  files  were  uploaded  to  a  Sun 
conversion  to  ASCII  tne  aaud 

sparcstation  for  analysis.  The  proper  time  base  was  added  to 
the  data  files,  and  then  the  oscilloscope  traces  contained 
within  the  data  files  were  transformed  into  freguency  domain 
via  a  discrete  Fourier  transform  so  that  corrections  for  cable 
and  probe  attenuation  could  be  made.  The  oscilloscope  trace 
data  was  then  discrete  inverse  Fourier  transformed  back  into 
the  time  domain,  plotted,  and  examined  for  data  quality.  Data 
from  shots  in  which  the  pulse  power  device  did  not  perform 
properly,  and  data  in  which  the  oscilloscopes  did  not  trigger 
properly,  were  eliminated  from  the  data  set. 

Circuit  modelling  software  had  been  used  to  estimate  the 
rising  and  falling  voltage  slew  rates  to  be  expected  for  a 
circuit  of  this  nature.  Data  which  did  not  exhibit  falling 
voltage  slew  rates  in  the  neighborhood  of  these  estimates  were 
eliminated  from  the  data  set  since  a  falling  slew  rate  of 
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insufficient  magnitude  indicated  that  a  discharge  had  not 
occurred. 

The  problem  of  establishing  the  beginning  of  the  pulse  had  to 
be  addressed.  For  the  purposes  of  this  investigatron  the 
beginning  of  the  pulse  was  defined  to  occur  when  the 
electromagnetic  field  exceeds  the  do  breakdown  value  of  15.54 
V/cm/Torr.  This  definition  is  due  to  the  pulse  width  of  the 

qince  the  conditions  necessary 
applied  electromagnetic  field,  since  one 

for  dc  breakdown  are  never  achieved. 

The  theory  upon  which  this  investigation  is  based  assumes  that 
a  triangular  shaped  pulse  excites  the  nitrogen  to  breakdown 
(refer  to  appendix  B) .  Data  which  exhibited  waveforms  other 
than  triangular  in  nature  were  eliminated  from  the  data  set. 
Por  this  purpose  the  criteria  used  re<raired  that  no  localrzed 
minima  exist  between  the  beginning  of  the  pulse  and  the  pornt 
at  which  breakdown  occurred.  The  resulting  data  yield  afte 
elimination  of  suspect  shots  was  approximately  48.6  percent. 

TWO  Channels  of  optical  data  were  acguired  to  provide  a 
verification  of  the  N,  spectrum  produced  by  arcing.  A  cable 
25  m  long  viewed  the  discharge  perpendicular  to  the  arc  while 
a  cable  50  m  long  viewed  the  discharge  down  the  axis  of  the 
arc.  The  optical  data  was  acguired  by  the  E.E.aG.  O.M.A.  and 
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stored  to  data  files  in  a  format  peculiar  to  the  O.M.A.  The 
data  files  were  converted  to  ASCII  by  a  software  program 
developed  at  the  Phillips  Laboratory,  and  then  uploaded  to  a 
Sun  SparcStation.  The  data  was  adjusted  for  fiber  optic  cable 
attenuation  and  geometrical  effects  of  probe  placement,  and 
then  plotted  and  examined  for  the  characteristics  of  the  N, 
spectrum.  The  data  was  not  adjusted  for  CCD  sensitivity  since 
the  effects  of  the  sensitivity  were  egual  in  both  optical  data 


channels . 
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Ill  D.  DATA  ANALYSIS 


The  data  which  was  acquired  as  described  in  the  experimental 
procedures  section  is  presented  in  figure  21  on  page  50,  a 
plot  of  E„/P  versus  Ptp.  The  ordinate  represents  E./P  in 
V/cm/Torr  while  the  abscissa  represents  Pt,  in  Torr-s.  The 
solid  line  is  the  theoretical  relativistic  estimate  of  Graham 
and  Roussel-Duprd,  while  the  dashed  line  is  the  classical 
estimate  as  discussed  in  the  classical  theory  section.  The 
symbols  represent  the  breakdown  points  acquired  during  the 
course  of  experiment  which  survived  the  data  reduction 

process . 

The  data  from  each  pressure  setting  except  for  0.336  Torr  to 
0.790  Torr,  and  0.930  Torr  was  averaged  to  obtain  mean  values 
and  standard  deviations  for  E„/P  and  PX,.  Only  the  PXp  values 
obtained  for  0.336  Torr  to  0.790  Torr,  and  0.930  Torr  were 
averaged.  The  E„/P  values  obtained  for  0.336  Torr  to  0.790 
Torr,  and  0.930  Torr  were  not  averaged  as  the  data  at  this 
pressure  establishes  the  lower  limit  of  PXp  under  the 
prevailing  experimental  conditions . 
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50 


Classical  and  Relativistic  Theory 


The  mean 


values  of  E./P  and  PXp  with  their  standard  deviations 
Shown  in  brackets  are  presented  below  in  table  2  for  each 

.operational  test  pressure. 


pr-occ..  (Torr)  Mean^ZP  (y/cm/Tor^ 


685 

295 

150 

73 


29.5 

19.5 
7.0 
4.8 
3.4 


2.08 


1.22 


1.07 

1.01 

0.890 

0.790 

0.700 

0.590 

0.506 

0.506 

0.336 


0.337x10^ 

0.638x10^ 

0.112x10^ 

0.160x10^ 

0.165x10" 

0.162x10" 

0.304x10^ 

0.623x10^ 

0.923x10^ 

0.135x10" 

0.201x10" 

0.228x10" 

0.189x10" 

0.165x10" 


(0.060x10^) 

(0. 056x10’) 

(0.007x10") 

(0.011x10") 

(0.018x10") 

(0.037x10") 

(0.028x10") 

(0.100x10’) 

(0.313x10’) 

(0.045x10") 

(0.072x10") 

(0.063x10") 

(0.072x10") 

(0.039x10") 


(18kV) 

(20.125kV) 


0.930  (Transformer  Rebuild) 


Table  2.  Mean  E^/P 


PT,.  (Torr-^ 


0.145x10'^ 

0.706x10-' 

0.421x10-' 

0.216x10-' 

0.794x10-’ 

0.589x10-’ 

0.211x10-’ 

0.134x10-’ 

0.983x10-® 

0.516x10-® 

0.286x10-® 

0.239x10-® 

0.200x10-® 

0.207x10-® 

0.200x10-® 

0.176x10-® 

0.150x10-® 

0.119x10-® 

0.118x10-® 

0.162x10-® 

0.102x10-® 


(0.030x10-®) 
(0.048x10-') 
(0.009x10-') 
(0.003x10-') 
(0.065x10-’) 
(0 . 046x10"’) 
(0.023x10"’) 
(0.025x10"’) 

(0.086x10"®) 

(0.018x10-®) 

(0.040x10"®) 

(0.012x10-®) 

(0.040x10"®) 

(0.023x10'®) 

(0.007x10-®) 

(0.010x10-®) 

(0.008x10-®) 

(0.008x10-®) 

(0.002x10-®) 

(0.012x10-®) 

(0.006x10-®) 


and  PXp. 


standard  deviations  of  the 
Data  which  did  not  lie  within  ±  3  stanaa 

„ean  values  was  regarded  as  suspect  and  elimnated  from  the 
population.  The  data  which  regained  is  presented  in  figure  22 
on  page  52,  a  plot  of  E„/P  versus  Pt.  The  ordinate  represents 
E  /P  in  V/cm/Torr,  the  abscissa  represents  Px,  in  Torr-s.  The 
s°olid  line  represents  the  estimate  of  Graham  and  Roussel- 
Dupr4,  while  the  dashed  line  is  the  classical  estimate. 
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Figure  22.  Edited  E„/P  vs.  PXp  Data 
compared  to  Classical  and  Relativistic  Theory 


V,  an  values  with  error  bars  representing  ±1.5  standard 
The  mean  values,  wx  •  .  n  n-sn 

.d-.  ^  Torr  to  0.790  Torr,  and  0.93U 

deviations,  and  the  edited  0.336  Torr 

in  fiaure  23  on  page  54,  a  plot  of 
Torr  data  are  presented  m  g 

versus  PV  -  ordinate  represents  E„/P  in  V/c./Torr 

while  the  abscissa  represents  Pt,  in  Torr-s.  Asm 

hhe  solid  line  is  the  relativrstic  estimate  while  the  das  e^ 
line  is  the  classical  estimate.  The  symbols  represen 
„ean  breakdown  points  or  edited  data,  calculated  as  desorr  e 

above . 

.he  discontinuity  in  the  region  ol  approxrmately  3xl0-»  Torr-s 

be  due  to  multichanneling  ot 
to  4x10-®  Torr-s  is  suspected  to  be  du 

.  Tt  has  been  observed  in  the  MFX  experiment 

the  discharge.  It  has  oeen  u 

that  at  approximately  30  to  40  Torr  arcing  phenomeno 
transitions  from  single  channel  arcing  through  multichanne 
arcing  to  a  symmetric  aural  discharge,  and  then  transitions 
hack  to  single  channel  arcing  at  approximately  3  to  4  Torr  ^ 
The  WFX  experiment  was  an  investigation  of  the  utilisation 
hydrogen  gas  as  a  working  fluid  for  the  production  of  plasma^ 
with  the  multichanneling  effect  behaving  as  a  smoo 


transition  from  40  Torr  to  4  Torr. 
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Figure  23.  Mean  and  Edited  E./P  vs.  Pt, 
compared  to  Classical  and  Relativistic 


The  multichanneling  effect  would  be  physically  limited  by  the 
nurrtoer  channels  available.  The  arc  diameter  of 

experiment  as  measured  from  Polaroid  photographs  of  single 
arcing  events  indicates  a  maximum  of  0.5  cm  per  channel.  The 
width  of  electrode  available  for  arcing  is  7.23  cm,  which 
allows  approximately  14  simultaneous  channels.  Assuming  that 
each  Channel  has  the  same  resistance  yields  an  eguivalent 
resistance  for  an  aura  discharge  which  is  one  fourteenth  that 

of  the  single  channel  resistance.  This  is  illustrated  in 

T  ^  t-Tna)-  \T  ths  inult ichsnnBlinQ 

figure  24  (A)  and  (B)  .  Note  that  V„,  the 
voltage,  is  the  square  root  of  one- fourteenth  V. 


Figure  24.  Multichanneling  Geometries 
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For  a  one  ns  RF  pulse  this  would  indicate  a  transition  region 
on  the  paschen  curve  from  approximately  3x10-’  Torr-s  to 
approximately  4x10-*  Torr-s.  A  transition  region  such  as  this 
is  supported  by  the  experimental  data,  which  shows  a  drop  in 
the  E„/P  values  in  the  transition  region  by  roughly  one  order 

of  magnitude. 


The  data  was  corrected  for  the  multichanneling  effect  by 
graphically  selecting  an  E,/P  value,  and  the  corresponding  Pt„ 
value,  of  greatest  magnitude  above  4xl0->  Torr-s  and 
graphically  selecting  the  most  comparable  E„/P  value,  and  the 
corresponding  Pt,  value,  below  4xl0-‘  Torr-s.  The  lowest  value 
of  E./P,  and  the  corresponding  Pi,  value,  in  between  these  two 
boundary  values  was  assumed  to  approximately  be  the  location 


of  the  symmetrical  aural  discharge.  These  assumptions  were 
made  due  to  the  uncertainities  introduced  because  of  the 
difference  in  the  working  fluids  of  each  experiment,  hydrogen 
and  nitrogen.  The  Pt,  values  selected  for  the  upper  and  lower 
bounds  were  7.94x10-*  Torr-s  and  2.39x10-’  Torr-s  respectively. 


The  symmetrical  aura  was  assumed  to  occur  at  2.11x10  Torr  s. 
A  straight  line  was  visually  fit  to  the  data  lying  between  the 
aural  discharge  and  the  lower  bound  via  a  log-log  graph.  The 


data  was  then  adjust  by  a  logarithmic  amount  to  compensate  for 
the  multiohannelling  effect.  Vl8  was  selected  as  the  maximum 
scaling  factor  since  a  symmetrical  aura  could  support  a 
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minimum  of  the  equivalent  of  18  simultaneous  channels. 
Likewise,  a  straight  line  was  visually  fit  to  the  data  lying 
between  the  aural  discharge  and  the  upper  bound  via-  a  log-log 
graph.  The  data  was  then  adjust  by  a  logarithmic  amount  to 
compensate  for  the  multi channel ling  effect  again  scaling  the 
result  to  yield  a  value  of  VlS  at  2 . 11x10  Torr— s . 

The  Eq/P  values  and  their  standard  deviations  were  rescaled 
yielding  a  Paschen  curve  as  shown  in  figure  25  on  page  59,  a 
plot  of  Eo/P  versus  PTp.  The  ordinate  represents  E^/P  in 
V/cm/Torr  while  the  abscissa  represents  PXp  in  Torr-s.  The 
solid  line  is  the  theoretical  relativistic  estimate  of  Graham 
and  Roussel -Dupre,  while  the  dashed  line  is  the  classical 
estimate  as  discussed  in  the  classical  theory  section.  The 
symbols  represent  the  breakdown  points  corrected  for 
multichanneling . 

Having  corrected  for  the  multichannel  effects  the  data 
conforms  most  precisely  to  the  classical  theory  developed  in 
this  investigation,  with  the  exception  of  a  shift  in  the 
magnitude  of  E„/P.  This  shift  in  magnitude  is  indicative  of 
an  incorrect  value  for  Vejf.  Estimation  of  an  appropriate  Veff 
yields  a  value  which  is  approximately  6  times  that  classically 
assxjmed.  The  classical  theory  developed  in  this  investigation 
was  adjusted  to  account  for  the  greater  value  for  v^ff  and 
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recomputed.  The  results  of  this  reevaluation  are  shown  in 
figure  26  on  page  60,  a  plot  of  Eq/P  versus  Ptp.  The  ordinate 
represents  E<,/P  in  V/cm/Torr  while  the  abscissa  represents  PXp 
in  Torr~s .  The  solid  line  is  the  theoretical  relativistic 
estimate  of  Graham  and  Roussel-Dupre,  the  long  dashed  line  is 
the  classical  estimate  as  discussed  in  the  classical  theory 
section,  while  the  short  dashed  line  is  the  classical  estimate 
corrected  The  symbols  represent  the  breakdown  points 
corrected  for  multi channeling. 


in 

o 


(a:ioi/uio/A)  a/03 
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Figure  25.  Corrected  Eq/P  vs.  PXp  Values 
Compared  to  Classical  and  Relativistic  Theory 


A  tabulation  of  the  mean  pulse  length  to  breakdown  at  each 
pressure  is  presented  in  table  3 .  Also  presented  is  the 
number  of  points  which  constitute  the  mean.  The  overall  mean 
pulse  length  to  breakdown  was  2.373  ns. 


Pressure  in  Torr 

0.930 

685 

295 

150 

73 

29.5 

19.5 
7.0 
4.8 
3.4 
2.08 
1.22 
1.07 
1.01 
0.890 
0.790 
0.700 
0.590 
0.506 
0.506 
0.336 

Table  3  .  Bre 


•  in  ns 

#  of  Pts 

1.103 

11 

2.117 

10 

2.393 

10 

2.807 

10 

2.959 

10 

2.692 

10 

3.201 

7 

3.014 

7 

2.792 

5 

2.891 

18 

2.481 

9 

2.344 

6 

2.234 

5 

1.980 

6 

2.326 

6 

2.532 

3 

2.513 

10 

2.551 

8 

2.304 

5 

2.332 

7 

1.638 

24 

own  Pulse 

Length 

Comparing  the  experimentally  obtained  spectra  acquired  by  the 
spectrograph  during  a  discharge  to  the  spectra  for  nitrogen 
from  the  CRC  Handbook  suggests  that  the  light  pulse  produced 
by  the  electrical  discharge  is  due  to  electron  energy 
transitions  in  nitrogen.  No  aluminum  lines  were  evident  in 
the  spectra. 
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A  discrete  Fourier  transform  was  perform  on  the  data  to 
compensate  for  the  attenuation  of  the  cables  and  probes.  A 
typical  Fourier  transform  of  an  electrical  discharge  is  shown 
in  figure  27  on  page  63,  a  plot  of  frequency  versus  amplitude 
after  compensation  for  cable  and  probe  attenuation.  The 
ordinate  represents  frequency  in  Hz  while  the  abscissa 
represents  amplitude  in  V/Hz.  The  localized  maxima  in  figure 
27  represent  the  maximum  probe  attenuation  at  those 
frequencies,  refer  to  figure  17.  The  overall  peak  frequency 
component  -  neglecting  the  localized  maxima  -  occurs  at 
approximately  100  to  200  MHz. 

The  previous  efforts  utilized  predominately  monofrequency 
sources.  The  results  of  Felsenthal  and  Proud  were  obtained 
for  a  pulsed  D.C.  experimental  system,  and  the  results  of  Ali 
and  Coffey  were  calculated  for  10  GHz  and  3  GHz.  Thus  the 
broadband  signal  should  be  bounded  on  the  bottom  by  the  work 
of  Felsenthal  and  Proud,  and  on  the  top  by  the  results  of  Ali 
and  Coffey.  This  is  graphically  presented  in  figure  28  on 
page  64,  a  plot  of  E^/P  versus  PTp.  The  ordinate  represents 
Eq/P  in  V/cm/Torr  while  the  abscissa  represents  PTp  in  Torr-s. 
The  boundary  limits  apply  between  approximately  4x10'®  Torr-s 
to  5x10'^  Torr-s,  where  departures  from  the  limiting  values 
occur.  These  departures  are  attributed  to  the  broadband 
nature  of  this  investigation. 
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Discrete  Fourier  Transform  of  N2  Discharge 
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Figure  28.  Corrected  Eq/P  vs.  PXp 
Corrected  Classical  Theory  Compared  to  Prior  Investigations 


CONCLUSIONS 


The  data  suggests  that  multichanneling  occurs  between  2.39x10'^ 
Torr-s  and  7.94x10'®  Torr-s.  This  implies  that  gaseous  switch 
design  may  be  modified  to  take  advantage  of  the 
multichanneling  effect  to  produce  switches  which  close  at  a 
lower  Eq/P  value  by  allowing  for  multiple  channels  or 
symmetric  aura.  This  further  implies  that  gaseous  switch 
design  may  be  modified  to  take  advantage  of  the 
multi channeling  effect  to  produce  switches  which  can  close  at 
selectable  E^/P  values  for  a  set  PipValue. 

The  data  corrected  for  multichanneling  suggests  that  the  v^ff 
assiimed  for  classical  monofrequency  Paschen  curve 
investigations  underestimates  the  for  broadband  Paschen 
curve  investigations  by  a  factor  of  approximately  6.  This 
would  imply  that  the  is  frequency  dependent. 

The  reversal  feature  of  the  Paschen  curve  has  not  been 
observed  in  the  experimental  results.  The  relativistic 
Paschen  curve  of  Graham  and  Roussel-Dupre  may  not  have  been 
noted  below  4x10'®  Torr-s  for  the  following  reasons: 

The  relativistic  theory  has  been  developed  to 

describe  phenomena  in  the  upper  atmosphere.  This 
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implies  a  theory  based  on  fixed  pressure  and 
variable  pulse  length.  The  experimental  data 
acquired  in  this  investigation  is  based  on  variable 
pressure  and  fixed  pulse  length.  It  has  been 
proposed^^  that  the  divergence  from  the  relativistic 
theory  to  the  classical  theory  at  approximately 
4xl0‘®  Torr-s  could  be  a  transition  region  from 
fixed  pressure  and  variable  pulse  length  physics  to 
variable  pressure  and  fixed  pulse  length  physics. 

The  relativistic  theory  of  Graham  and  Roussel-Dupre 
naturally  requires  that  the  electrons  achieve 
relativistic  velocities.  This  issue  been  examined 
in  the  context  of  this  investigation^^.  A  velocity 
estimate  had  been  obtained  by  dividing  the 
acceleration  of  a  free  electron  in  a  1x10’  V/m  field 
by  the  classical  estimate  for  at  1  Torr, 
yielding  a  P  of  approximately  1.  To  calculate  the 
collisional  velocity  of  this  investigation  requires 
a  knowledge  of  the  mean  free  path  length’^. 
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1 


(25) 


X- 


N„Ttk  ij 


where;  =  debye  length 

Nn2  =  Density  of  nitrogen 
Referring  to  equation  C-25  it  is  recalled  that 

~  ^  (26) 
23 .8x10^Pq\  0 


Using  the  classical  physics  equation  x  =  ^t^, 
recalling  the  Lorentz  force  equation  niea  =  qE^jeid  to 
calculate  acceleration,  and  setting  x  equal  to  the 
mean  free  path  length  yields  an  equation  for  the 
mean  collisional  time,  tcoii*  Numerically  solving 
for  tcoii/  and  then  calculating  the  electron  velocity 
at  the  time  of  collision  yields  an  estimated 
maximum  electron  velocity  of  P  =  1.7x10'®.  This  is 
insufficient  for  the  observation  of  relativistic 
effects . 

The  theory  of  Graham  and  Roussel-Dupre  has  been 
developed  assuming  intense  E-fields,  B-fields,  and 
propagation  vectors.  The  experimental  setup  of 
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this  investigation  reflected  an  RF  pulse  off  of  an 
open  transmission  line,  effectively  canceling  the 
incident  B-field  and  propagation  vector  with  the 
reflected  B-field  and  propagation  vector,  see 
figure  29.  The  relativistic  fluid  eguations^^ 
collapse  to  classical  equations  under  these 
conditions.  For  relativistic  effects  to  be 
observed  requires  either  a  B-field  or  electron  flow 
in  the  direction  of  the  propagation  vector. 


Tip  Sensor 


Figure  29.  Field  Geometries  at  Tip 

Classical  theory  and  experimental  data  suggest  that  a 
collisional  frequency  limit  of  approximately  lxl0‘®  Torr-s 
exists.  This  implies  that  an  intense  electrical  pulse  may  be 
created  in  a  gaseous  medium  without  breakdown  occurring  if  the 
parameter  of  pulse  length  multiplied  by  pressure  does  not 
exceed  1x10'®  Torr-s. 
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APPENDIX  A 


Evaluation  of  the  Ionization  Rate 
Integrals 
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Evaluation  of  integral  11  as  listed  in  the 


main  text 


400 


|^e-V2[a  J  E(E-E^)exp{^^}exp{^}S  + 

P  2e 


102  f  VEexpi-^jdE] 

400 

may  be  performed  in  the  following  manner. 


<~>-6xl0 


27 

2z 


400 


f-3/2 


[aexp{-^}  f  E^exp{-[^^-^]EidE~ 

Ei  ^ 


400 


as^expl-^}  f  Eexp{-  EfdE 


10 


^  f  y/Eexp{-^}dE] 

Anfx  « C 


A-2 


Evaluate  each  integral  separately.  The  strategy  to  be  used 
will  begin  with  the  evaluation  of  the  second  integral  in 
equation  A-2,  continue  with  the  first  integral,  and  conclude 
with  the  evaluation  of  the  last  integral. 
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First  evaluate  the  second  integral  in  equation  A-2 .  Begin  by 
defining  k„  =  1/p  +  3/2e  =  [2e+3p] / (2Pe) . 

400  400 

f  Eexp{-k^£!tdE-^  J  exp  i-k^)  dE  A-3 

El  O 


_ 3  expC-ic^g)  ^00  g  r  exp  (-400k J -exp  (-F.kJ 

K  ‘‘  ‘ - r, - 


Thus 


400 

J  Eexp  ( -kgE)  dE —  [  (-400jCo)  -exp  (-E^kg) 

Eh  kl 


^  -400e3q)  (-400jCo)  +E^exp(-E^k^)  , 


A-5 


However  «  400  eV  and  1/k,  =  2Pe/[2e+3P],  therefore 
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A- 6 


400 

/ Eexpi-k^E)dE-.^^^SzEi^A.  [A+£-.] 


and 


400 


}  Eexp(-[±*-^]^dE- 


23e 

2e+3p 


zHi 

2e 


) 


A-7 
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Evaluation  of  the  first  integral  in  equation  A-2  is  best 
accomplished  by  taking  the  partial  with  respect  to  of  the 
second  integral  in  equation  A-2.  Thus,  the  partial  with 
respect  to  k^  of  equation  A-6  will  yield  the  solution. 

400 

S^exp  { -k^)  dE-  f  Eexp  ( -k^E)  dE^  o>  [  1  i 

ok^  i  °  dk^  k^  k^  ^ 


/ 


^2exp{-E^k^)  ^  2E^6xp{-E^k^)  ^  Ejexpi-E^kg) 


kl 


Thus 


“expi-E^kg)  [-^  + 

■Kq 


] 


400 

J  E^exp  ( -kgE)  dE^ 


expC-gjkp) 

K 


+e!] 


and  finally 


400 

J  E^exp{-  [  -g- + ]  j^dg“ 


2pe 

2e+3P 


r  8(Pe)^ 

(2E+3p)2 


4  peg. 
^  2e+3P 


+g|]  exp( 


)exp( 


-3gi 

2e 


) 


A-8 
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Finally  evaluate  the  third  integral  in  equation  A-2,  with  the 
definition  of  k  =  3/2e,  and  the  substitution  of  x  = 

M  m 

xBxp ( -kx^ ) dx^ 

400  20 


-2  fx^expl-kx^)  dx--2-J-  r exp(-kx^)  dx~-2^l 
20 


Now  note  that 


••  2k 


fexpl-k(x^+y^)]cixdy-  f  f rexp(-kr^)  drdB 


2020 


■ysoTT  0 


•• 

J  rexp{-kr^)dr-Tz  J  exp{-kr^)dr^ 


(-icr2)^oo-^exp  (-800ic) 


■N 


exp  (-400k) 

JC 


74 


Thus 


j*  /ffexp{ 


J- 1-2  3 


400 


2e 


dk 


dk\ 


exp  (-400ic) 

jZ 


J  y/Eexp{ 


1 


ejqj  (-400ic)  +800a 


-^exp(-400ic) 


ultimately  yielding 


j  v^exp{-^^}d£:-* 

400  ' 


27ce 

3 


exp{-600/e)  [-^+800] 


A- 9 
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The  evaluation  of  the  integrals  in  equation  A-2  therefore 
become 


<  —  >-6xl0\ 

P  >271 


2aPe 

2e+3P 


r  8(pE)2  ,  4pegj 

(2e+3P)2  2e+3p 


+£'|]exp( 


-3gi 

2e 


) 


2aPe£'^ 

2e+3p 


_2£|_.B^)exp(C2^) 


+10‘ 


2^exp(^^)  (21*800)) 
3  e  3 


A-10 
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Appendix  B 


Derivation  of  the  Effective  E-Field 
for  Short  Electromagnetic  Pulses 
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Consider  the  effective  electromagnetic  field  for  continuous 
wave  phenomenon® . 


^eff” 


_ E _ 

<ff 


B-1 


It  is  desired  to  obtain  the  effective  electromagnetic  field 
for  a  short  electromagnetic  pulse,  which  contains  a  spectrum 
of  frequencies.  For  simplicity  assume  that  the 
electromagnetic  pulse  is  a  sawtooth  waveform  as  shown  in 
Figure  B-1. 


Figure  B-1.  E-field  vs.  Time 
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This  waveform  is  described  by 


where : 


for 


E(t) 

Eq  =  Peak  electric  field 
t  =  time 

tp  =  pulse  duration 


B-2 


and  zero  elsewhere.  The  Fourier  series  of  this  sawtooth 
waveform  is 


“  9  c>  . 

^sin{-52Lt)  B-3 

n-l  tp 

where  njc/tp  is  the  component  frequency,  and  ( -1 ) '’■"^2Eo/n7t  is  the 
E-field  due  to  each  frequency  component.  The  effective  E- 
field  at  each  frequency  component  is  obtained  by  substitution 
of  the  E-field  due  to  each  frequency  component  of  equation  B-3 
into  B-1.  This  operation  yields. 


‘'P 


im 


2^0 

[2  (l  +  -il!lEL)  ]  1/2 
4^lff 


B-4 
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To  obtain  the  total  effective  E-field  it  is  necessary  to  sum 
over  all  frequencies. 


iJ-l 


7371  (1  +  -^!^)  1/2 

x-2  2 

off 


B-5 


Recalling  once  more^°  that  Vejj=5 . 3xlO®P  yields 


i2-l 


_ y/2  (-1)”^! _ 

nn  (1+ - - )i/2 

2.809xl0^^tpP^ 


B-6 


The  summation  of  the  first  100  terms  contained  in  equation  B-6 
is  plotted  in  Figure  B-2,  on  page  81,  as  a  function  of  pulse 
length  multiplied  by  pressure.  The  asymptotical  value  for  the 
summation  is  approximately  0.310  for  the  region  of  interest  of 
this  investigation.  Therefore,  equation  B-6  may  be  rewritten 
as 


E^f^-0.310E^ 


B-7 


At  higher  pressures  this  value  is  achieved  for  shorter  pulse 
lengths,  while  at  lower  pressures  a  longer  pulse  length  is 
required.  Figure  B-2  indicates  that:  a  large  applied  E-field 
can  be  scaled  into  an  appropriately  small  effective  E-field 
with  the  selection  of  an  adequately  short  pulse,  and  as 
pressure  decreases  so  does  the  effective  E-field. 
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aniBA  uoT^Bunnns 


81 


Figure  B-2 .  Summation  Value  vs.  Pulse  Length 


Appendix  C 


Relativistic  Derivation  of  the 
Paschen  Curve 

for  Short  Electromagnetic  Pulses 
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Consider  relativistic  scattering  by  a  central  force,  as  shovm 
in  figure  C-one  below. 


Figure  C-One.  Scattering  by  a  Central  Force, 
where:  b  =  impact  parameter 

dQ  =  differential  solid  angle  for  a  unit 
radius  sphere  (  27crsin0d0,  with  r=l  ) 

FinS  incident  particle  flux 

Tout-  particle  flux  scattered  into  dQ  at  0 

da  =  particle  flux  scattered  into  dQ  at  0,  divided 

by 

Following  the  method  of  Cherrington^^ ,  the  differential  cross 
section  can  be  defined  as' 


do  particle- flux- into-dQ/dCl  i  ^^outi 
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Conservation  of  particles  requires  that 


{2%bdb)  (27isin0d6) 


Therefore  via  substitution, 

da  _  ^ouci  _  2‘Kbdh _ b  |  cib| 

dQ  r,„  ®  27rsin6d6  sin0  d0 


The  relationship  between  b  and  0  is 


71-0 

2 


-/ 


du 


21  1/2 


C-2 


where:  u  =  1/R 

E  =  kinetic  energy 
V  =  potential  energy 

The  electric  field  for  a  relativistic  particle  is 

E _ 0 _ 

47ce 

where:  Q  =  electron  charge  (  1. 602x10’^®  Coul . ) 

Eq  =  electric  permittivity 
P  =  v/c 

Y=(l-p2)-i/2_ 
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The  potential  energy  for  a  repulsive  interaction  then  becomes 


V-qr<p  (i2,4))- 


_ _ 

47te  (1-P^sin^<|)) 


C-3 


Noting  from  figure  C-1  that  sinO=b/R,  equation  C-3  may  be 
rewritten  as 


V- 


_ o! _ 

47ce^Y2[l-(-^)2]"/" 

iv 


For  compactness  of  notation  the  quantity  k  is  now  defined  as 


k~ 


_o!_ 

Arte^ 


so  that 


_ k _ _ _ uic _ 

y2jJ[l_(_&fc)2]3/2  y^[l- 

R 


Equation  C-2  may  now  be  written  as 


7t-6 

2 


_ du _ 

{  1  Uk{v) _ 

2j2  Y^E(l-{pijU)2)3/2 


]  -u2}1/2 


C-4 
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To  evaluate  the  integral  in  equation  C-4  it  is  necessary  to 
determine  u„ax-  This  is  accomplished  by  beginning  with  the  law 
of  conservation  of  momentum. 

which  may  be  rewritten  as 

yvb  _ _ 0 _ 

c 


Solving  for  the  time  derivative  of  (])  yields 


<P- 


yvb 


i?2 


N 


14(XE!*)2 

R 


C-5 


which  reduces  to  the  classical  answer  when  v«c,  such  that  p=0 
and  y=l .  The  law  of  conservation  of  energy  requires  that 


(Y-l)m^c2-(Yjj-l)/noC2+(Yjj^-l)m^c2+ - £_ - 

R^ 


C-6 


where 


Yj<- 


R  \  2 


l-(-) 

c 
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and 


1 


2 


C-7 


The  quantity  of  interest,  11^3^,  is  evaluated  at  that  point  in 
space  at  which  the  radial  velocity  is  zero. 


This  drives  the  first  term  on  the  right  hand  side  of  equation 
C-6  to  zero.  Thus  equation  C-6  may  be  rewritten  as 


y2i?(l-l!^)3/2 


yzn^c^~yj^^m^c^+ 


Y2i?{l_l!^)3/2 


Combining  C-5  and  C-7  into  the  above  equation  yields 


R^ 


C-8 
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It  is  noted  that  the  classical  solution  for  goes  to  1/b  as 
velocities  increase.  This  implies  that-  a  relativistic 
solution  should  yield  approximately  the  same  answer.  To 
simplify  the  solution  of  equation  C-8,  the  second  term  on  the 
right  hand  side  is  approximated  for  R=b  and  P=1  as 

_ k _  k  _ky^  _yk 


Therefore  equation  C-8  becomes 


i+( 


jk 

R 


Solving  for  1/R  yields 


^^ax■ 


2y^k^ 

(2Y^iC)2  yk  )2_(Ypi3)2]  (yZ-l) 

2[(-li^)2-(YPi3)2] 

-i-±  (— (Y^-l)  [(— 

\  Y^ _ _  C-10 

[(-^)"-(PJb)2] 
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Noting  that 


A(Y2-1)-P2 

y2 


Equation  C-10  becomes 


^nax* 


(_^^)2-p2[{^iL_)2-(P2))2] 


m^c‘ 


[( 


)"-(P2))2] 


which  can  be  rewritten  as 


^nax“ 


(■ 


[( 


ittgC^pb 


)2-l] 


The  denominator  is  functionally  equal  to  minus  one. 

( - - - )^+  — 

yntgC^fi^b^  b^ 

which  yields  the  classical  solution  for  v«c,  and 
u™ax=l/b  for  v=c. 


^nax*’ 


nigC^^^b^ 


Thus 

C-11 

yields 
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It  is  also  possible  to  obtain  the  critical  impact  parameter, 
b^^it/  and  hence  calculate  the  inelastic  collisional  cross 
section,  b^^it  is  defined  to  be  the  solution  for  b  when  is 
negative  or  complex.  Recalling  the  equation  prior  to  C-11 


[( 


)2]-l 


^nax 


C-12 


( 


— - )2+  — 


which  yields  a  b^rit  of 


^crit^ 


P- 


p/n^c 


C-13 


Thus  the  inelastic  cross  section  becomes 


o(P) 


inBlastic“'^^cTi  t 


Tlk^ 


(/noC2)2p2 


C-14 


Rewriting  in  terms  of  energy 


O  inelastic-T^b. 


2 

'crlt’ 


71  iC^ 


Kc2)2[i- 


1  +  - 


E 


C-15 
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Uniax  having  been  obtained,  the  integral  C-4  may  now  be 
evaluated. 


7t-6 

2 


^DdX 


-/ 


du 


_ _ 

Y2E(l-(piDU)2)3/2 


]  -u2}l/2 


Recall  from  equation  C-9 

_ k  ■  k 


therefore 


71-6 

2 


du 


{—i - U-U  2)1/2 

i?2  E2)2 


C-16 


This  equation  is  merely  the  classical  equation  for  a  repulsive 
coulomb  interaction  with  the  second  term  of  the  denominator 
modified  by  the  relativistic  y.  Immediately  utilizing  the 
classical  results,  and  replacing  the  classical  "k”  with  the 
relativistic  "yk"  yields  the  solution 


7C-6 

2 


--sin  ^ 


[-2U- 


yk 

Eb^ 


] 


yk  \2]  1/2 
Eb^ 


C-17 
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The  term  contained  by  the  first  arcsin  on  the  right  hand  side 
of  C-18  goes  to  -1  as  v-^c.  This  allows  C-18  to  be  reduced  to 


—  -sin' 
2 


jrt) 

\  2-]  1/2 

Eb2 


which  can  be  written  more  compactly  as 

0-2tan-i-^  C-19 

Equation  C-19  reduces  to  the  classically  obtained  result  when 
y— >1 .  Solving  for  b  yields. 

Recalling  C-1 

do  b  I  dt»|  1  ,  Y-*^\2  1 

dQ  “  sine  dO  ”  4  '  2E  ^  sinM  — ) 
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The  integral  for  the  collisional  cross  section  then  becomes 


eisinM|) 


C-21 


which  yields 


.(e)-rt(^)"cot2(5|ii) 


C-22 


where  0„in  is  0  evaluated  at  b  set  equal  to  the  species  debye 
length 

The  solutions  for  and  Gg  are  also  desired,  where  G^,  is  the 
forward  momentum  transfer  cross  section,  and  g^  is  the  energy 
cross  section.  These  parameters  are  defined  as  followed. 

[1-COS  (0)  ]  27isin{0)  d0 

®Biln 


o_(e)-  f  -^[i-cos(e)]  (^)27isiii(e)d0 

J  aU  M 

®mtn 
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The  classical  solutions  immediately  show  that 


o„(E)  -1671  ( 


yk 

4E 


)nn( 


2EXd 

Yic 


C-23 


o^(E) 


2^ 

M 


o^(E) 


C-24 


Thus,  the  solutions  for  all  three  collision  cross  sections 
contain  the  species  debye  length.  The  species  debye  length  is 
itself  subject  to  relativistic  effects.  The  debye  length  is 
defined  to  be 


\  1/2 

Q 


where:  n^,  =  electron  density  in  particles  per  cubic  meter 
kg  =  Boltzmann's  constant 
Te  =  electron  temperature 
kgTe  =  electron  kinetic  energy  in  eV 

The  relativistic  kinetic  energy  is 

E-Cy-D/WoC^ 
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Thus , 


At  breakdown  the  ratio  of  the  final  electron  density  to  the 
initial  electron  density  has  been  determined^’  to  be  10®.  At 
approximately  standard  temperature  and  pressure  the  free 
electron  density  has  been  estimated’®  to  be  between  0.1  to  10 
electrons  per  cubic  meter.  This  investigation  will  assiame  a 
free  electron  density  of  one  electron  per  cubic  meter . 
Utilizing  the  ideal  gas  law,  PV=n„oieRoT  where:  P  is  pressure  in 
atmospheres,  V  is  volume  in  liters,  n^oie  is  the  number  of 
moles,  Rq  is  the  gas  law  constant  of  0.08205  atm*l/mole/°K,  and 
T  is  the  temperature  in  degrees  Kelvin,  the  initial  electron 
density  n^  may  be  estimated  as 

®  V  RoT 

where  the  factor  of  23.8x10®  in  the  numerator  is  in  units  of 
electrons*liters  per  mole*ms’.  Substituting  in  for  n^ 

C-25 


23. BxlO^PQ 
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Substituting  equation  C-25  into  . equation  C-19  yields  a  of 


e„,.„-2tan 


-1  Yic 

2E^ 


23 .SxlO^Q^P 


Therefore  equation  C-22  evaluates  as 


23. SxlO^Q^P 


C-26 


which,  except  for  the  energy  term,  is  independent  of  y. 
Equations  C-23  and  C-24  become 


) 


C-27 


®  Af  4E  yk\  23. BxlO^Q^P 


C-28 


Noting  that  y=  l+E/m„c^  C-27  maybe  rewritten  as 


o,(E)-7i[^(l  +  -^)]nn(4J - 


k\  23  .QxlO^PQ^  1 1  E  C-29 


m^c‘ 


while  equation  C-28  is  still 
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Therefore,  with  the  use  of  equation  7  of  the  main  text  and  a 
standard  text  reference^^. 


V 

where  C-31,  C-32  and  C-33  all  have  units  of  s'^.  Substitution 
of  C-15,  C-26,  and  C-29  into  C-31,  C-32  and  C-33  respectively 
yields 


J.  r^exp(^) 

-2e 


C-33 


V  , -6 .02x1026^ 


27ir  ,  k  X2-T 


)H-l 


3-  Eexp(-3^) 


2e 


in,  - 1  ] 


-dE 


1+- 


E 


m^c‘ 


C-34 


6 .02x1026 


23.8xl0®^| 


277t  (££)e’2|E2exp(_^) 


dE 


C-35 


V  -6 .02x102® 


2711 


C-36 


/-i(l+^)nn( 

J  il/  m.c^ 


k\ 


- ^  . -^-  ^  ^ -- )  exp  ( ^ )  dE} 

23.8X10®P0^  1 ,  E  2e 


in,c' 


Evaluating  integral  C-34  requires  numerical  integration. 
Dividing  the  integral  by  pressure  in  Torr,  and  equating  with 
equation  5  of  the  main  text  yields 
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Evaluating  integral  C-35  yields . 
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Integral  C-36  requires  numerical  integration.  From  equation 
17  of  the  main  text  it  is  recalled  that 


C-39 


It  is  desired  to  obtain  the  relativistic  effective 
electromagnetic  field  for  a  short  electromagnetic  pulse,  which 
contains  a  spectrum  of  frequencies.  Let  the  applied  pulse 
shape  be  that  selected  in  appendix  B  and  begin  by  recalling 
equation  B-5  which  is  still  correct  in  the  relativistic 
regimen . 
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It  has  been  shown^^  to  be  inappropriate  to  assume  that 
Veff=5 .3xlO®P  with  the  ratio  of  the  applied  electric  field  to 
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pressure  is  in  excess  of  50  V/cm/Torr.  Roussel-Dupre  et .  al.^’ 
has  noted  that  for  intense  electric  fields,  which  produce 
relativistic  motion  in  electrons,  the  effective  collision 
frequency  is 

where  R.  =  10*n„v./V,„-  Rc  =  10“n„v„/V,„,  and  V,„  is  the  volume 
of  the  arc. 

Equations  C-36  through  C-41  constitute  a  set  of  existential 
equations.  Since  the  experimental  parameters  of  E^,  and  P  are 
typically  known,  the  following  strategy  for  the  estimation  of 
the  breakdown  curve  was  adopted.  Veff  was  first  set  equal  to 
the  classical  value  of  5.3xlO®P,  tp  was  set  equal  to  1x10'^  s, 
and  Eq  was  set  equal  to  200  V/cm.  These  values  were  then  used 
to  estimate  E^fj  based  on  equation  C-40.  From  this  value  for 
Eeff,  an  estimate  of  e  was  obtain  in  accordance  with  equation 
C-39.  Utilizing  this  estimate  of  e,  an  estimate  of  Vi  and 
hence  tp  was  obtained  via  equation  C-37.  If  the  estimate  of 
the  pulse  length  differed  from  the  previous  estimate  of  the 
pulse  length  by  greater  than  one  percent,  then  a  new  estimate 
of  Egff  was  obtained  and  the  process  was  iterated.  If  the 
value  for  the  pulse  length  had  converged,  then  estimates  for 
Vc  and  Vro  were  obtained  via  equations  C-3  6  and  C-38.  A  new 
estimate  for  was  then  calculated  via  equation  C-41.  If 
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the  new  estimate  of  V^ff  differed  from  the  previous  estimate  by 
greater  than  one  percent,  then  a  new  estimate  of  was 

obtained  and  the  entire  process  was  iterated  until  convergence 
to  a  final  value  was  observed. 

The  results  of  this  computation  are  shown  in  figure  C-2,  a 
plot  of  Eq/P  versus  PXp.  The  ordinate  represents  E^/P  while 
the  abscissa  represents  PTp.  The  solid  line  is  the 
■  theoretical  curve  of  Graham  and  Roussel-Dupre^^ .  The  long 
dashed  line  is  the  theoretical  curve  as  derived  in  the 
classical  theory  section  of  the  main  body.  The  short  dashed 
line  is  the  theoretical  curve  produced  by  the  process 
described  in  this  section. 

Note  that  all  three  theories  produce  a  reversal  in  the  paschen 
curve.  In  the  classically  derived  theory  the  reversal 
transitions  into  a  vertical  barrier.  In  the  relativistic 
theory  of  Graham  and  Roussel-Dupre  the  reversal  transitions 
through  yet  another  reversal  to  reestablish  the  Paschen  curve 
at  a  higher  value  of  E^/P.  The  relativistic  theory  developed 
in  this  section  produces  a  reversal  which  continues  into  the 
upper  right  hand  quadrant  of  the  plot,  which  is  attributed  to 
the  one  dimensional  nature  of  the  theory. 
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